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We measured Shubnikov-de Haas (SdH) oscillations in GaAs (100) two-dimensional holes to deter- 
mine the inversion asymmetry-induced spin splitting. The Fourier spectrum of the SdH oscillations 
contains two peaks, at frequencies /_ and /+ , that correspond to the hole densities of the two spin 
subbands and a peak, at frequency /tot, corresponding to the total hole density. In addition, the 
spectrum exhibits an anomalous peak at /tot/2. We also determined the effective masses of the two 
spin subbands by finding the inverse transform of the Fourier spectrum in the vicinity of /_ and 
/+, and then analyzing the temperature dependence of the SdH oscillations for each subband. We 
discuss our results in light of self-consistent calculations and previous experiments. 

PACS numbers: Valid PACS appear here 



In a solid that lacks inversion symmetry, the spin-orbit 
(SO) interaction leads to a lifting of the spin degeneracy 
of the energy bands, even in the absence of an applied 
magnetic field. In such a solid, the energy bands at finite 
wave vectors are split into two spin subbands with dif- 
ferent energy surfaces, populations, and effective masses. 
This so-called zero-field spin splitting is particularly im- 
portant for semiconductor holes because they typically 
have much smaller kinetic energies than electrons. The 
problem of inversion asymmetry-induced spin splitting in 
two-dimensional (2D) carrier systems in semiconductor 
hetero junctions and quantum wells Q, |^ |^ has become 
of renewed interest recently 3 because of the possible 
use of such systems in realizing spintronic devices such 
as a spin field-effect transistor 0, Q, and for studying 
fundamental phenomena such as the spin Berry phase 
L?, 8]. 

GaAs 2D holes in modulation-doped GaAs/AlGaAs 
heterostructures 0, Q provide an attractive system for 
these studies as they can have high low-temperature mo- 
bilities and a strong SO interaction that is tunable via 
the application of gate bias P, T^l . Recent work in these 
systems has concentrated primarily on 2D holes that are 
fabricated, via molecular beam epitaxy (MBE), on GaAs 
(311)A substrates. The main advantage of using (311)A 
substrates, rather than the more commonly used (100) 
substrates, is that one can use Si as a p-type dopant in 
the MBE growth. The standard p-type dopant for MBE 
growth on GaAs (100) is Be. Unlike Si, however. Be has 
a tendency to diffuse and migrate under typical MBE 
growth conditions and, as a result, fabrication of high- 
quality GaAs (100) 2D holes has been challenging. On 
the other hand, because of a higher symmetry, the band 
structure for (100) GaAs 2D holes is simpler than for 
(311)A j^. It is therefore desirable to study the prob- 
lem of spin splitting in GaAs (100) 2D hole systems. We 
have been able to fabricate such 2D systems with reason- 
ably high mobilities, and report here our measurements 
of their spin splitting and of the hole effective mass, and 



the comparison of the experimental data with the results 
of state-of-the-art energy band calculations |3, ^3 • 

The holes in GaAs (100) were in fact the first 2D 
system in which the inversion asymmetry-induced spin 
splittingwas observed experimentally using SdH oscil- 
lations 0, . The oscillations in 2D holes confined to 
GaAs/AlGaAs heterojunctions revealed a beating pat- 
tern, indicating the presence of two spin subbands with 
different hole densities 2\ . The spin splitting was inferred 
from a change in the frequency of the oscillations above 
a certain magnetic field. The temperature dependence of 
the amplitude of the oscillations along with the assump- 
tion that the spin subbands are parabolic was used to 
deduce their masses 0|- Although the results of these 
experiments were in qualitative agreement with the cal- 
culations that followed, some questions remain, particu- 
larly regarding the effective mass values jl^i, ilSj iMi llJI ■ 
In our study, we revisit the spin splitting in GaAs (100) 
2D holes via a careful analysis of the beating pattern 
in SdH oscillations. We employ Fourier transform (FT) 
techniques to determine the temperature dependence of 
the amplitude of the SdH oscillations and the effective 
masses for both the lighter (HHl) and heavier (HHh) 
heavy-hole spin subbands independently. Our method 
has the advantage that it does not require assumptions 
regarding the band structure, namely, the parabolicity of 
the spin subbands for determining their effective masses. 
The masses we find are in good agreement with the re- 
sults of subband calculations. The data and calculations 
also provide clear evidence for the strong nonparabolicity 
of the HHh band. 

Our samples were grown on GaAs (100) substrates 
by MBE and contain modulation-doped 2D hole sys- 
tems that are confined to either a 20 nm-widc GaAs 
square quantum well or a GaAs/AlGaAs heterojunc- 
tion. The square well is flanked on each side by undoped 
Alo.3Gao.7As spacer layers. On the front (surface) side, 
a Be-doped layer (Be concentration of 2.6 x 10^* cm~'^) 
of Alo.sGao.r follows the 21 nm thick Alo.sGao.? spacer 
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FIG. 1: (a) Observed Shubnikov-de Haas oscillations for a 2D hole system confined to a (100) GaAs/AlGaAs heterojunction. 
(b) The Fourier spectrum of the oscillations (in the range O.ST < B > 3T) at two different temperatures. The dashed boxes 
show the square windows applied to isolate the /_ and /+ peaks before obtaining the inverse Fourier transforms which are 
shown in the insets in Fig. |21 (c) Fourier spectrum of calculated magneto-oscillations. The open circles indicate the expected 
peak positions according to the calculated spin subband densities at zero magnetic field (see text). 



layer. The Alo.3Gao.7As/GaAs interface in the hetero- 
junction sample is separated from a 16 nm thick Be- 
doped Alo.3Gao.7As layer (Be concentration of 3.5 x 
10^* cm~^) by a 25 nm Alo.3Gao.7As spacer layer. In or- 
der to reduce Be diffusion and migration during the MBE 
growth, and hence increase the quality of the samples, the 
substrate temperature was lowered to 550°C from 640°C 
before doping. In addition, for the square well sample, 
the doping was done on the front side of the well only 
unlike previous samples, which had doping on both the 
front and the substrate sides. We fabricated Hall bar pat- 
terns on all the samples and used In/Zn alloyed at 440°C 
for the ohmic contacts. Metal gates were deposited on 
the front and the back of the samples to control the den- 
sity. The typical low temperature mobility for the square 
well sample is 2.6 x 10^ cm^/Vs at a 2D hole density of 
p = 1.7 X 10^^ cm^^ and for the heterojunction sample it 
is 7.7 X 10"* cm^/Vs at p = 2.3 x 10" cm'^ Longitudi- 
nal (Rxx) and transverse (Rxy) magneto-resistances were 
measured, as a function of the perpendicular magnetic 
field (B), at T « 30 mK via a standard low frequency 
lock-in technique. 

Figure^a) shows the low-field SdH oscillations for the 
heterojunction sample. The beating pattern in the oscil- 
lations indicates that the two spin subbands are nonde- 
generate. The FT spectrum of the oscillations, shown 
in Fig. n^b), exhibits four dominant peaks at frequen- 
cies /avg, /+, and /tot, with the relation /tot = 
/+ + /- = 2 /avg. The /tot frequency, when multiplied 
by e/h, matches well the total 2D hole density deduced 
from the Hall resistance (e is the electron charge and h 
is Planck's constant). The two peaks at /_ and /+ cor- 
respond to the SdH oscillations of the holes in individual 
spin subbands although, as discussed below, their posi- 
tions times e/h do not exactly give the spin subband den- 
sities. As we also discuss later in the text, the presence 
of a peak at /avg is related to this anomalous behavior. 



We note that the /avg peak is observed in all our GaAs 
(100) 2D hole samples, including the square well sam- 
ple which shows an FT spectrum qualitatively similar to 
Fig. ^b) . In the remainder of the paper we concentrate 
on the data from the heterojunction sample as it shows 
a much higher spin splitting than the square well. 

We first compare the measured FT spectrum of 
the SdH oscillations with that of calculated magneto- 
oscillations of the density of states (DOS) at the Fermi 
energy (Ep) We obtain the Landau fan chart at 

B > by evaluating an 8 x 8 k • p Hamiltonian that fully 
takes into account SO coupling due to both the struc- 
ture inversion asymmetry of the GaAs/AlGaAs hetero- 
junction as well as the bulk inversion asymmetry of the 
underlying zinc blende structure [Til IT^ . In our calcula- 
tions we assumed that the system formed an accumula- 
tion layer. For such systems Stern pointed out that 
the Hartree potential depends on the concentration of mi- 
nority impurities. We assumed that the concentration of 
unintentional minority impurities in the GaAs layer was 
1 x lO^"' cm~'^ and we used a spacer width of 25 nm; both 
assumptions are consistent with our sample parameters. 
We note that one should not expect a one-to-one cor- 
respondence between the measured SdH oscillations and 
the calculated magneto-oscillations of the DOS at Ep- 
For example, the amplitude of the oscillations which de- 
termine the amplitudes of the FT spectrum peaks are 
different from the SdH oscillations. However, the peak 
positions in the FT spectrum - the quantities we are in- 
terested in here- are not affected by these details [T^ITt) . 

Figure ^c) shows the FT spectrum of the calculated 
DOS magneto-oscillations for the sample parameters of 
Figs, n^a) and (b), i.e., a heterojunction sample with 
p — 2.3 X 10^"'^ cm^^. The spectrum shows three main 
peaks, marked by /_, /+, and /tot- The positions of 
these peaks are in good agreement with the correspond- 
ing three peaks in the experimental data [Fig. ^b)]. 
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The /avg peak is also seen, albeit with a smaller am- 
plitude, in the FT spectrum of the calculated magneto- 
oscillations. In Fig. njc) we have indicated, by open cir- 
cles, the positions of FT peaks expected from the relation 
/ = {h/e)p±, where p± are the densities of the two spin 
subband densities calculated at B = 3- It is evident 
that the frequencies /_ and /+ do not directly give the 
spin subband densities via the relation p± = {e/h)f±. 
This anomalous behavior, like the unexpected peak at 
/avg J reflects the fact that the spin precession along the 
cyclotron orbits becomes nonadiabatic if the SO coupling 
in some parts of ky space is sufficiently weak Quali- 
tatively, one can argue that the peak at /avg is due to the 
presence of cyclotron orbits in which the holes move half 
the time along the constant energy contour of the HHh 
band and then switch to the HHl contour (16] . This sim- 
plified picture cannot explain, however, that the peaks at 
f± do not obey the relation f± — {h/e)p±. We remark 
that anomalous magneto-oscillations were observed also 
in GaAs (311) A 2D electrons and holes [Tll20| . 

We now concentrate on our results for the effective 
masses of the HHl and HHh spin subbands. In a 2D 
system with only one spin subband occupied, or with de- 
generate subbands, the variation in the amplitude AR^x 
of the SdH oscillations with temperature T is commonly 
used to determine the carrier effective mass m*. This 
is done by fitting AR^x to the Dingle factor, £/ sinh£, 
where ^ = 27^^ ksT / {hujc) and = eB/m* HHH. Here 
771* is the fitting parameter. In our system however, due 
to the spin splitting of the subbands giving rise to beating 
patterns in the SdH oscillations, the two masses cannot 
be deduced directly by fitting ARxx to a Dingle factor. 
In the SdH experiments performed by Eisenstein et al. 
|2| on GaAs (100) 2D holes, the authors assigned the 
low field oscillations (with frequency /_) to the SdH ef- 
fect of the spin subband with the lighter mass, and the 
oscillations at higher fields (with frequency /tot) to the 
total 2D hole density. They then determined the mass for 
HHl, 7771, by fitting the temperature dependence of the 
low field oscillations to a Dingle factor. To deduce the 
HHh mass, 777^, they assumed that the two spin subbands 
have parabolic dispersion curves i?±(k||) = ti^ k'^^ / {2m'^) , 
where ky is the in-plane wave vector and m^. are con- 
stants independent of the energy E. This model implies 
the relation m'^/m*^ = /-/(/tot — ./-) which was used in 
Ref. to obtain m*_^_. 

In our study, the FT spectrum of the oscillations allows 
us to independently determine rrij. for the two spin sub- 
bands. We isolated the individual peaks, /_ and /+ for 
HHl and HHh, with square windows [Fig. ^b)] and cal- 
culated their inverse FT The insets in Figs. [2ta) 
and (b) show the inverse transform for each of these 
peaks. The SdH oscillations in these insets can then be 
attributed to the HHl and HHh spin subbands, respec- 
tively. The main parts of Fig. |2] show the temperature 
variation of the amplitude of these oscillations. We fitted 
this variation at particular B values to a Dingle factor for 
each subband [Figs. Ela) and (b)] and deduced the 777* 
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FIG. 2: Variation of ARxx with temperature and its fit to the 
Dingle factor with m* as the fitting parameter for (a) the HHl 
and (b) the HHh subbands, respectively. Insets in (a) and (b) : 
After isolating the peaks /_ and /+ by square windows [see 
FigH^b)], inverse FT is performed on these peaks separately 
to deduce SdH oscillations of the HHl and HHh subbands. 



values as the fitting parameters. Figure O shows 777* val- 
ues determined as a function of B for a hole density of 
2.3 X 10^^ cm^^. The amplitude variation of the oscil- 
lations fitted to the Dingle factor at different values of 
B should in principle give the same mass. In our sys- 
tem however, ttiI has a strong, nearly linear dependence 
on the value of B at which the variation of Rxx with 
temperature is analyzed. Only 777 is approximately in- 
dependent of B. The origin of the field dependence of 
m*_ which was also seen in Ref. remains unknown. 
Linear extrapolation of the data to B = (dashed lines 



in Fig. 13), suggests 



0.9 7770 and m* 



0.2 7770, 



where mo is the mass of free electrons. Note that the ra- 
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FIG. 3; Values for ml and m!j. for the HHl and HHh sub- 
bands, deduced from fits of SdH oscillation amplitudes to 
the Dingle factor for the heterojunction sample, with p = 
2.3 X lO^'^ cm~'^ (mo is the free electron mass). Inset: m* 
calculated from the subband structure at -B = 0. 



tio of the hole densities of the two subbands is not equal 
to the ratio of their masses. Hence, in general it cannot 
be assumed that the bands are parabolic and, in partic- 
ular, m*f_ cannot be deduced from the ratio of the spin 
subband densities and ml. 

We determined m* at different hole densities ranging 
from 2.0 to 2.4 x 10^^ cm~^ for the heterojunction sample, 
where we used front and back gate biases to vary the 
density. The qualitative trends were similar in all cases. 
The strong dependence of m*_ on B was observed at all 
densities and was found to be relatively independent 
of B. Quantitatively, m"^ ranged from 0.75 to 0.91 and 
ml extrapolated to a _B = value between 0.18 to 0.25 
in this density range. 

We compare the measured effective masses with an 
energy-dependent DOS effective mass 



ml{E) 1 



mo 



n 2mo 



(1) 



for the two spin subbands. Equation |^ gives the effec- 
tive mass m*j^{E) one needs locally at energy E to ap- 



proximate i?±(k||) by a parabolic dispersion which yields 
the same DOS like E +Ck n). Our numerical calculations 
of i;±(k||) follow Ref. [nj. We then evaluate Eq. Q by 
means of analytical quadratic Brillouin zone integration 
|2^ . We remark that the SdH oscillations are an effect 
taking place at Ep- Therefore, the DOS m*, defined by 
Eq. is the appropriate quantity for comparison with 
the effective mass deduced experimentally from the tem- 
perature dependence of the SdH oscillations. However, 
one should not expect to find a one-to-one correspon- 
dence between experiment and theory because the tem- 
perature dependence of magneto-oscillations in 2D sys- 
tems can be more complicated than implicitly assumed 
by the Dingle formula [25ll2fil|. 

The inset in Fig. O shows the DOS effective mass cal- 
culated for p = 2.3 X 10^^ cm~^. The subband edge is 
at —17.9 meV and Ep at —19.4 meV. We obtain an ef- 
fective mass at Ep of 0.9 toq and 0.2 toq for the HHh 
and HHl bands, respectively. This is in good agreement 
with our experimental results (Fig. (SJ. From the curve 
for the heavier mass, we see that the HHh band is highly 
nonparabolic. We note that symmetry requires that in a 
Taylor expansion of the subband dispersion i?±(k||), the 
-6 = spin splitting is characterized by terms propor- 
tional to odd powers of ky . Only the spin-independent 
terms in the Taylor expansion are proportional to even 
powers of ky [23. In 2D hole systems the odd-power 
terms give rise to a spin splitting that is comparable in 
magnitude to the spin- independent part of the subband 
dispersion |13. Ila ll4L Il5 | . The spin-split subbands thus 
have nonparabolic dispersion curves with strongly spin- 
dependent masses. We remark that spin-split electron 
subbands can usually be characterized by the same value 
of m* because in the electron case, spin splitting rep- 
resents only a small correction to the spin-independent 
part of the subband dispersion • 

In conclusion, we have measured SdH oscillations in 
GaAs (100) 2D holes. Through Fourier analysis of the 
oscillations, we determined the splitting as well as the 
effective masses of the two spin subbands. The results 
are in good agreement with the calculations. There re- 
main however, some open questions regarding the origin 
of the average peak observed in the FT spectrum and the 
dependence of experimentally determined mass values on 
B, especially for m*_. 
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